Auger decay of degenerate and Bose-condensed excitons in Cu 2 



ON 



X 



G. M. Kavoulakis and Gordon Baym 
Department of Physics, University of Illinois at Urbana- Champaign, 1110 West Green Steet, Urbana, Illinois 61801 

(Submitted to Physical Review B, 15 February 1996) 

We study the non-radiative Auger decay of excitons in CU2O, in which two excitons scatter to an 
excited electron and hole. The exciton decay rate for the direct and the phonon-assisted processes is 
calculated from first principles; incorporating the band structure of the material leads to a relatively 
shorter lifetime of the triplet state ortho excitons. We compare our results with the Auger decay 
rate extracted from data on highly degenerate triplet excitons and Bose-condensed singlet excitons 
in Cu 2 0. 



I. INTRODUCTION 



Recent experiments on excitons in CU2O have been carried out at sufficiently high densities and low enough 
temperatures that the Bose-Einstein statistics they obey becomes significant and, indeed, in stressed samples, 
Bose-Einstein condensation of excitons has been observed. As we have shown in Ref. 0], the crucial barrier to 
condensing excitons is heating of the exciton gas by the non-radiative Auger decay process, in which two excitons 
collide to produce an ionized electron and hole. See Fig. 1. This process not only heats the gas, but it leads to loss of 
excitons as well. [The Auger process is familiar in electron-hole plasmas, in which the recombination of an electron 
and a hole excites either an electron high in the conduction band or a hole deep in the valence band. This process has 
been studied there both theoretically || and experimentally Q.] In this paper we present detailed calculations of the 
Auger process in CU2O that we used in Ref. Q to calculate the time evolution of the exciton gas in CU2O following 
laser excitation. 



In the CU2O experiments, intense pulses of laser light excite the crystal, creating a gas of (triplet) ortho excitons 
and (singlet) para excitons, split by an exchange energy AE » 12 meV ||. The kinetic energy distribution of ortho 
excitons as a function of time from the onset of the laser pulse is observed by spectroscopy of their photoluminescence 
(considerably more intense than that of the para excitons). In the classical (low density) regime, the energy distribution 
is observed to be Maxwcll-Boltzmann, described by an effective exciton temperature. In the quantum (high density) 
regime, the spectra are well fit in terms of an ideal Bose-Einstein gas with an instantaneous chemical potential, fi, 
and temperature, T. From these parameters, one can directly calculate the density of the gas. The experiments of 
Snoke et al. observe a "quantum saturation" of the ortho excitons, i.e., a tendency for them to move closely parallel 
to the critical line without condensing The critical line, an adiabat (constant entropy per particle, s = S/N), has 
the form, T c = (2nh 2 /mk B )[n/g((3/2)} 2 / 3 , where the de generacy 5 = 3 for ortho excitons and 1 for para excitons. 
In general lines in the phase diagram parallel (in a log-log plot) to the condensation line at higher temperature are 
adiabats. 

The important role of the Auger process has been revealed in several experiments. The inverse of the total decay time 
of ortho excitons in CU2O shows an approximate proportionality to their density; the data for non-stressed crystals 
[^|,^) is shown in Fig. 2. In this figure the quotes indicate that the density plotted on the horizontal axis is assumed 
to be proportional to the total radiative recombination luminescence intensity and is calibrated by comparison to the 
density deduced from spectral fits. The total decay rate is a sum of interconversion of ortho excitons into para excitons 
H and the Auger process; the rates of both processes are comparable. An additional feature of the Auger process 
is that it leads to a much shorter lifetime of ortho excitons; in lightly stressed crystals the observed ortho-exciton 
decay time is ~ 0.1 ns Q compared with the para-exciton lifetime at high para-exciton density, ~ 10 19 cm' 3 , which 
exceeds 100 ns. The Auger process also leads to the presence of ortho excitons in quantum wells, caused by high 
stress, for times much longer than the ortho-exciton lifetime; although the exciton gas in the wells consists primarily 



1 



of para excitons, high energy electrons and holes produced in Auger ionization of para excitons lead to reformation 
of excitons in essentially random internal angular momentum states JToj ] . 

Here we calculate the Auger decay rate for both the direct and phonon-assisted mechanisms, Figs. 1 and 3. One of 
our basic conclusions is that the rate of the phonon-assisted Auger mechanism is much larger than the rate of the direct 
process, a consequence of the conduction and the valence bands having the same parity. Although the phonon-assisted 
mechanism requires participation of a phonon, the rate is enhanced because the process is dipole allowed. In Sec. 
II, we first derive the matrix element of the direct Auger process, neglecting for simplicity at this stage the angular 
momentum structure of the excitons, then we discuss detailed approximations to the matrix element and calculate the 
decay rate for the direct process. In Sec. Ill we generalize the results to take into account the full angular momentum 
structure of the ortho and para excitons in CU2O and explain how specific properties of the excitons in this material 
are related to the direct Auger process. We then turn, in Sec. IV, to consider the phonon-assisted Auger process and 
discuss its connection with the phonon-assisted radiative recombination of excitons in CU2O. In next section, V, we 
take into account effects of the Pauli principle on the decay rate. Finally, in Sec. VI, we compare our results with the 
measured rate and summarize our conclusions. 



II. DIRECT PROCESS - CALCULATION OF THE MATRIX ELEMENT AND THE DECAY RATE 

The direct Auger process is described by the sum of the topologically distinct Feynman diagrams shown in Fig. 1. 
In this figure, time progresses from left to right; the solid lines with forward-going arrows denote electrons, and those 
with backward-going arrows denote holes. The initial state contains two excitons, of center-of-mass momenta K and 
P, described by the wave function 

*i = * K (Xe,X h )*p(x^X^), (1) 

where denotes the full wave function of an exciton of total momentum Q, and the Xj denote electron and hole 
coordinates. We include effects of symmetrization below. The final state contains an ionized electron of momentum 
k e and an ionized hole of momentum kh, described by the wave function 

*/ = $c,k e (x e )$.,k h (x h ), (2) 

where the ^^(x) are the normalized Bloch wave functions for the ionized electron and hole in the final state, and the 
subscripts c, v denote the conduction and valence bands. The wiggly line in the figure denotes the screened Coulomb 
interaction. In the first class of diagrams, A and A', the electron and hole in the same exciton recombine; in the 
second class, B and B', particles of different excitons recombine. 

The exciton wave functions ^Q(x e , Xh) can be expressed in terms of the electron and hole Bloch wave functions as 

* Q (x e , X h ) = ^2<f>q * C ,q+Q/2(Xe) < I , t, -q+Q/ 2 (Xh), (3) 

q 

where the Bloch wave functions ^^(x) in band j are 

^,k(x) = ^, k (x)e*- x . (4) 

The momentum q in the <&i q in Eq. (^) is restricted to be inside the first Brillouin zone. In general the sum on 
the right of this equation should include all the conduction and valence bands, but for simplicity we neglect all other 
bands except the lowest conduction and the highest valence, which is a rather good approximation for the case of 
CU2O. In the effective mass approximation for the exciton wave function, with equal electron and hole effective masses 
(which is approximately true for the bands we consider), 'J'q takes the form 

* Q (x e , x h ) = -L e »Q-(*e+x h )/2 0rci(Xe _ Xh)) (5) 
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where ft is the crystal volume, and (j> ro i is the relative electron-hole wave function. Comparing Eqs. (||) and (g) we 
see that q is the Fourier transform of (j) re \ times VT2; to a good approximation, we may make this identification even 
when using complete Bloch waves, Eq. (^). 

The interaction between the excitons is, in general, the dynamically-screened Coulomb interaction, V Q s(k.,Lu) — 
Vbouiomb(k)/e(k, lu), where the dielectric function is given by the usual Lindhard result 

£ (k,c) = £oo - V(k) £ "l',k + q -"l.q | ( | >j9) (6) 

with the high frequency dielectric constant of the material, and I and I' band indices. The second term in (^) can 
in fact be neglected in the Auger process, where we are interested primarily in small momentum transfers of order 
the thermal momentum of excitons, P 2 /2M CXC on the scale of 1-10 meV, but rather large frequencies corresponding 
to energy transfer on the order of the energy gap, Uu> w E g — 2.17 eV. In this limit, the terms in Eq. (||) involving 
the same band vanish for finite to. In the same limit, for different bands, I ^ V, the last term on the right side of Eq. 
(||) gives a fc-independent contribution which can be ignored for low densities. Thus e(k, to g ) ~ too, and the effective 
interaction, essentially independent of u> for large u>, becomes 

47TP 2 

We note that as a consequence of CU2O being almost non-polar pT| , f M (w 6.46) and the static dielectric constant 
eo (~ 7.11) differ only slightly. 

Assembling these results, we find that the matrix element M = (^/l^l^i) between the initial and final states, for 
diagram A for example, is thus 

M A = ^2 ^-k h +P/2Kff(G-K)^ q (u„ )q _ K /2|u C) q + K/2>G(Wc,kekc,ke+K>G'<^,+k h +G+G',K+P) ( 8 ) 



q,G,G' 



where 



(-"„, q _K/2|u c ,q+K/2)G = J dx<, q _ K/2 (x)w Ciq+K/2 (x)e lG - x (9) 

is a generalized overlap integral |L2"[ between the conduction and the valence bands. The momenta p in </> p are 
restricted to the first Brillouin zone. The x-integral is over a unit cell, and the summation is over all reciprocal lattice 
vectors G, G'. Energy and crystal momentum conservation restrict this summation considerably. Since K and P are 
determined by the thermal motion of excitons they are negligible compared to k e and kh, which are determined by 
the scale of the energy gap E g and the reciprocal lattice vectors. In this limit the condition for conservation of crystal 
momentum in Eq. (@) becomes, 

k e + k h + G + G' » 0. (10) 
In the same limit the energy conservation condition, 

fc2,2 t2,2 

l3L + l±h^ Eq , (11) 
2m e 2m/! 

puts an upper bound on the quanity |k e + kh|. Numerically, |k e + kh\ae < 5.5, where ag is the lattice constant. Since 
for the smallest nonzero reciprocal lattice vector Go, G$ai = 6.28, the conservation laws are satisfied only for 

G + G' = 0. (12) 

Using the fact that for a relative Is electron- hole state, <^_ q = q , the final expression for the sum of the matrix 
elements for processes A and A' in Fig. 1 is, 
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M A + M A > = ^2 ^=ff( G ~ K )</ ) q( W 1 >,q-K/2l« C ,q+K/2)G <Sk e +k h ,K+P 
q,G 

X (0k h -P/2(Wc : kJtic,k e +K)-G - 0k e -P/2(^,k h |ti«,k h +p)-G) ■ (13) 

The first term in parentheses in Eq. ( |l3| ) corresponds to process A and the second to process A', which enters with a 
change of sign because the Coulomb exchange couples to the positively charged hole instead of the electron. 

In the second class of diagrams, processes B and B' are equivalent to processes A and A' with either the two 
electrons or the two holes interchanged. For this reason these processes enter with a minus sign relative to A and A', 
and we find 

M B + M B > = - y^^q ^k e +k h ,K+p 

q,G 

x (0k h -p/2V r o ff(k e + G — q — K/2)(-u t , !q _ K/2 |M Cike _ K >-G(w c ,kel M c,q+K/2>G 

-</ , k e -K/2V c ff(k h + G -q-P/2)(ti 0)kh _p|u c>q _p/ 2 )_G(w«,kj,K,q+P/2)G) • (14) 

We now consider the effects of the band structure of CU2O on the generalized overlap integral, Eq. (|9|). For G = 0, 
(^) is the usual overlap integral between the conduction and the valence bands [^3|. We note that for K = 0, the 
wave vector of condensed excitons, the overlap integral for G = vanishes as a consequence of orthogonality of Bloch 
wave functions between different bands. 

The excitons in CU2O in the condensation experiments under study are formed from electrons in the lowest con- 
duction band (Tq) and holes in the highest valence band (the "yellow series"). The fact that these two bands 
have even parity O] implies that 

«3,k(-x) = Uj,-k(x). (15) 
This result is readily seen in the tight-binding approximation for the Bloch wave functions for which 

$ j , k (x) = ^$ a , i (x-R)e ikR , (16) 

R 

where the summation is over all crystal points and the <f> a j(x) are the atomic orbitals which form the band with 
index j; for positive parity, 

* A k(-x) = <&i,-k(x), (17) 
from which ( |l5| ) follows. Thus, for the G = term in the sum over q in the matrix elements Ma and Ma>, 

^ < / ) q('" t ;,q-K/2l u c,q+K/2)G=0 = ^ <^q (lVq+K/2 |«c,q-K/2) G=0 , (18) 

q q 

which implies that the sum is an even function of K and is thus proportional to K 2 (since it vanishes for K = 0). We 
use k • p perturbation theory [ pH) to calculate the overlap integral in the above matrix element. For example, 

I \~l \ l Tl \ " ( M »,ql K / 2 'PK.q) | , nn , 
K,q+K/ 2 ) » K,q) + ~ ~ ~ ~ Kq>) ( 19 ) 

where m is the bare electron mass. The sum is over all the negative parity bands of CU2O. Of the ten valence and four 
conduction bands in this material, only the conduction band (c') that is ss 449 meV higher than the conduction 
band (c) and the one very deep valence band (v') ~ 5.6 eV below the valence band (v) have odd parity, Fig. 6. 
The contribution of the odd-parity conduction band is expected to dominate the sum because of the small energy 
denominator in Eq. (|l9|). Thus 

/ I X _ ( h \ 2 K>,ql - K / 2 • PKq) K.q|K/2 ■ pKq) , on , 

( "- K/2,W) ~ [m) £ (,, q - £ „ q )(e„ q -,, q ) • (20) 
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The G = term (which we denote by "normal" ) in the sum (|T^) for the matrix elements Ma + Ma> , in the limit 
K,P <C fc e j fcft, is equal to 

(Ma + MA')normaI = 77" — 77^7 (0k e -K-P/2 - 0k 6 -P/2) 0q (««,q-K/2|«c,q+K/2) <^k e +k h ,K+P (21) 
£oo q 

~ fS?2 K ' (^) £ <^ Kq-K/ 2 k,q + K/ 2 ) <5 ke+ k h ,K+P. 

Taking re i(a;) to be the hydrogenic ground state wave function with excitonic Bohr radius as (~ 7 A) and using Eq. 
(pp|), we estimate the above quantity as 

(M ii\>t\ ~ 2 7 7re 2 (k e a B ) ■ (Ka B ) ft 2 x ^ (K ■ p„,j)(K ■ p i|C ) 

+ M Abnormal 7j , 7] V2\3~ ~ 2 "F72? \7 r^k e +k h ,K+P, (2,2.) 

S2eoo (1 + (keaByy ui X 2 (e CjQ - £j,o)(£i,o - e«,o) 

where Pij is the matrix element of the momentum operator between the Bloch states of the bands i and j at the zone 
center. The matrix elements pjj can be extracted from experiment; if we assume, for example, that they are all of 
equal magnitude, then Ref. [|l6| implies that \pij\/ti » 0.13 A -1 . 

If K and P are ignored compared to k e and kj l7 in the case that G = 0, the right side of Eq. (|||) vanishes. The 
contribution of Umklapp processes to Eq. (H3h is, 



(M A + Af^) Umk l a pp W 2J Kff(G) </>q</>k e ("-u,q|u C ,q)G ( (u C ,k e |Ue : k e ) -G - ("o.kh |^K,k h ) G ) ^k e +k h ,K+P, (23) 
q,G#0 

which is nonzero, but very small. In the nearly free electron model, for example, with u„,k(x) = l/\/7Q and u Ci k(x) = 
e °' x /-\/n7, where J7 C is the volume of the unit cell and Go one of the six smallest reciprocal lattice vectors, the right 
side of Eq. (|2^) vanishes. 

We turn now to the processes B and B'. In the limit K,P <C k e ,ky L Eq. ( |l4|) reduces to 

M B + M B > ~ 0q<5k e +k h ,K+P (0k h Kff(k e + G - q)(u c ,k e K,q)G(?VqK,ke)-G 

q,G 

- 0k. ^(kh + G q)(u„ ikh |u C:q )_ G (M 1 , !kh |u t , :q ) G ) . (24) 

For G 7^ 0, in the same nearly free electron approximation, the right side of Eq. (|4|) vanishes. Since the overlap 
integrals are less than unity in magnitude, an upper bound of the G = term is given by 

(M B + M S /) norma i < 20 ke ^ q Kff(k e - q)(u t ,, q |w c ,k e > ^k e +k h ,K+p- (25) 



q 



The overlap integral above is approximately 



^t),q|^c,k e 



H V V ^ - — ' Pc '' ^ - Pi,v (26) 
m ) ( £ c,o - e»,o)(e*,o - e«,o) 



where g ranges approximately between A: e — 1/a; and /c e + 1/a;. Then, 

(Mb + MB') no rmal < 41 — h 77 7T77 ( — J 77 l Pc ''IIP^I — I 7 <5k e +k h ,K+P , (27) 

SZeoo [1 + (k e a B Yr \ m J , (s Ct0 - e it o)(e ii0 - e v<0 ) 



where the numerical factor is the result of a dimensionless integral times constants. We show below using Eq. ( J2? 
that processes B and B' contribute negligibly to the Auger decay rate. Only A and A' contribute significantly to 
the direct Auger mechanism; these processes do not allow condensed excitons to Auger recombine, since Ma + Ma> 
vanishes for K = (Eq. (f|)). 
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We next proceed to calculate the decay rate Tk,p for the direct process in which two excitons of momentum K and 
P collide. From the golden rule 

Tk p = y E l M l 2 ( J - n cke)(l - n vMh ) S(E K + E P - e c , ke - e vMh ), (28) 

k e ,k h 

where M is the matrix element for the process. The exciton energies are 

E Q =E g -E b + ^-f-, (29) 

where M oxc is the exciton mass, Ef, is the binding energy, as 141 meV including the exchange energy AE « 12 meV 
which is nonzero for the ortho excitons only (see Sec. III). We may neglect E b and AE compared to the energy gap 
Eg. The electron and the hole energies are 

h 2 kl n 2 k 2 , s 

£ ^= E ° + 2^> £ ^ = ^ (30) 

and the rij ^s are the electron and hole occupation numbers in the final ionized states. Since these final states have 
very high energy, we may neglect the n^k, in Eq. (p8|), as well as the inverse process. 
The decay rate per unit volume for the direct Auger process, Ya^/Q., is given by 



^r A , rf = -i(^) =i^/ K / P r K ,P = -^, (3i) 



n Ad n\dtJ Ad si^" 4 "" 11 ' 1, T A ,d 

where /k is the distribution function of excitons, n the density of excitons and TA,d the scattering time for the direct 
Auger process. 

Using Eq. ( p2] ) we estimate the contribution to the decay rate due to the (dominant) processes A and A' 

r- A \ « 0.77 x 2 17 . >f U1 )3 1B ^ f — VE ( 1 N!^ *V k * T »W»|, (32) 

where E g = h 2 k g /2fi, fi is the reduced mass for the electron and the hole, n cxc is the density of noncondensed (excited) 
excitons and T is the temperature of the exciton gas. We note that TA,d is inversely proportional to the density of 
noncondensed excitons. For the proportionality constant for excitons in CU2O, 



' A,d 



5 x 10 _4 n oxc Tns _1 , (33) 



with T measured in kelvin and n cxc measured in units of 10 18 cm -3 . As we will see this is decay rate is small compared 
to the rate for the phonon-assisted process, but it can contribute to para-exciton para-exciton Auger collisions (which 
are, as we will see, forbidden in the zero-stress case) in highly stressed crystals. We explain this effect in the next 
section. 

Finally Eq. ( p7| ) gives an upper bound for i~ A 1 d . B+B ,, i.e., the contribution of processes B and B' to t^, 

, ^, 2x41 2 /xe 4 (k g a B ) ( h V ^ f |p e ,«l|p«,«l V „„ 3 ( o A] 

A,d;B+B' - ^ 2% Z ela {1 + (^^4 [ maB ) { _ £ . o) ( £ . Q _ £y Q ) ) ™B- W 

For the parameters of excitons in CU2O, Eq. ( |34| ) gives T A 1 d . g+B , < 4 x 10~ 4 nns _1 , with n measured in units of 
10 18 cm -3 , which is negligible compared with the contribution of the processes A and A' as well as the most important 
processes of the phonon-assisted Auger mechanism. 
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III. SPIN STATISTICS - BAND STRUCTURE AND OPTICAL PROPERTIES OF CU 2 



In this section we review the band structure underlying the properties of excitons in CU2O and their direct Auger 
recombination process. These properties infuence the direct as well as the phonon-assisted mechanisms. We discuss 
the effects of the band structure on the phonon-assisted Auger process in Sec. IV. 

The direct radiative recombination of both ortho and para excitons is dipole forbidden because the conduction and 
valence bands have the same parity jl4|,[l^] . The quadrupole recombination of ortho excitons is allowed and has been 
observed experimentally p8[ . In contrast, the direct radiative recombination of para excitons is highly forbidden [ jl9| , 
but becomes allowed if uniaxial stress is applied pc[ |. Both species can recombine via phonon-assisted processes. 

The exchange interaction, shown in Fig. 4, where an exciton virtually annihilates and reforms, is responsible for 
the ortho-para energy splitting at the zone center of CU2O. Although for pure spin states this interaction is nonzero 
for the singlet, shifting the energy of the singlet state higher than the triplet, experimentally the ortho excitons lie 
higher than the para excitons by 12 meV [EJ. The absorption spectrum of CU2O Jl7],|2l| shows discrete lines below 
the continuum which correspond to excitonic absorption, but the n=l line is very weak, in contrast with the n=2,3,... 
lines. All these observations can be understood in terms of the band structure of CU2O |l4| ]. The conduction band 
is formed by Cu 4s orbitals and the valence band by Cu 3d orbitals. The fivefold degenerate (without spin) Cu 3d 
orbitals split under the crystal field into a higher threefold and a lower twofold degenerate band r^ 2 . Finally, 
^25 splits further because of the spin-orbit interaction into two bands, a higher non-degenerate band and a lower 
twofold degenerate band (Fig. 5). The Ty" and Tg mix with the conduction band to form the yellow and green 
exciton series, respectively. 

From now on we restrict our discussion to the yellow exciton series. The weakness of the n=l line in the absorption 
spectrum is due to the Tq and bands having the same parity, making the transition to this line only quadrupole- 
allowed. In contrast to the n=l line, the n=2,3,... lines are dipole-active, provided that the created excitons are in 
p-states. The total angular momentum functions for the exciton triplet states are 



|J=1,J, = 1) =| Te.Tff) (35) 

\J=hJ z = 0) =^= (I Te, Ih) - I ie, Tff» (36) 
\J = 1,Jz = -1)=\L,Ih), (37) 



and for the singlet states, 



I J = 0, J z = 0) = i= (| te, in) + I |e, U)) • (38) 

The indices e, H refer to the electron and the hole, respectively: while the electron states are pure spin states, the 
hole states are total angular momentum states, 

\U) = -^l(X + iY)\l h )+Z\U)} (39) 

and 

\l H ) = -^[{X-iY)\W-Z\lh)], (40) 

where the states with lower case h are pure spin states. The spatial functions X, Y, Z transform as yz, xz and xy, 
respectively. 

Using these angular momentum functions we explain now the experimental observations mentioned above, turning 
to the direct Auger decay process at the end of the section. In the radiative recombination of an exciton the electron 
goes from the conduction to the valence band with emission of a photon. If we assume that there is no spin-flip (a 
higher-order effect), the above angular momentum functions imply that the matrix element for the direct radiative 
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recombination rate for the ortho excitons is proportional to v/2/3 times the result of the spatial part of the calculation, 
but that it vanishes for the para excitons. Similarly, para excitons are not allowed to recombine directly and the 
exchange interaction vanishes for the singlet. For the triplet state the exchange interaction is nonzero and raises the 
ortho-exciton energy with respect to the para-exciton energy; the energy splitting AE at the zone center (K = 0) is 
given by || 

AE=^J * K =o(x,x)^(x-x')*k=o(x',x')dxdx', (41) 

where V(x — x') is the Coulomb potential; the factor of 2/3 again comes from the angular momentum wave functions. 
Using Eq. (^) for the exciton wave function we find that 

AE =^ Y] ^q(^qK,qWpK,,pK,p)G. (42) 

The sum is dominated by the terms with smallest G's (six in number because of the cubic symmetry of the crystal), 
which we denote by Go; then 

A^<^^^|^. o1 (0)| 2 ^|Ko| Wc ,o)g o | 2 -^4— (— VlKok,o)Go| 2 . (43) 

Since experimentally AE « 12 meV [p2| , we have |(u„ i o|w c ,o)goI ~ 0.62. 

In Auger decay, processes B and B' are not affected by the angular momentum of the electron and the hole. The 
forbiddenness of the direct radiative recombination of the para excitons due to the spin conservation implies that the 
recombination vertices which appear in processes A and A' of Fig. 1 are nonzero only for an ortho exciton decaying 
and ionizing with either an ortho exciton or a para exciton. 

Uniaxial stress mixes the r^" and valence bands with the results that the para-exciton energy increases for 
small values of the applied stress but decreases for higher values, and the direct radiative recombination increases 
quadratically with increasing stress pp| . In the Auger process this mixing implies that the recombination vertex of 
stressed para excitons becomes allowed. Para excitons in highly-stressed crystals are therefore able to Auger recombine 
(processes A and A'), ionizing either an ortho exciton or a para exciton, just like the ortho excitons in the unstressed 
case. We discuss this effect further in Sec. VI. 



IV. PHONON-ASSISTED AUGER PROCESS 

The phonon-assisted Auger mechanism with the participation of a longitudinal, odd-parity optical phonon is in fact 
the dominant Auger process. Although this process requires a phonon and its rate is therefore reduced by factors 
of the exciton-phonon interaction, the dipole matrix element between the intermediate state and the conduction or 
the valence band does not vanish, in contrast to the case of the direct Auger mechanism. Since the matrix element 
for normal processes (G = 0) does not vanish for K — 0, the contribution of Umklapp processes should be relatively 
small, and we neglect them. In this case the angular momentum of the colliding excitons does not put any restriction 
on Auger collisions. Following the same procedure as before, we calculate all Feynman diagrams (Fig. 3). We assume 
for simplicity that the lattice temperature is very low, and only spontaneous phonon emission is possible. 

For the processes C, C and D, D' of Fig. 3 with two excitons of momentum K and P colliding, giving an electron 
of momentum k e , a hole of momentum kh and a longitudinal optical (LO) phonon of momentum Q and energy hujQ, 
denoted by the dashed line, the matrix element is 

M c + M c > + M D + Mjy = U off (K - Q) (0 kh _ P/2 - ke _ P/2 ) 

!Mt),q-K/2| u n,q+K/2-Q) (u n ,q+K./2-Q \HhO |w C ,q+K/2) 

>k e +k h +Q,K+Pj (44) 



£ c,q+K/2 — £«,q+K/2-Q — Tu*)Q 
'j,q-K/2|-ffLo|lVq-K/2+Q)(u„ iq ~K/2+Q|Uc,q+K/2) 



TlLO, 



Q 
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where i?LO is the carrier (i.e., electron or hole) - LO-phonon interaction. The first term in the sum refers to processes 
C and C, and the second to processes D and D'. The sum over n includes the two odd-parity bands of CU2O. We 
again use k • p perturbation theory to calculate the overlap integrals in the above matrix element. For example, 

|(Q-K). Pl u ri.q-K/2/ 

(«« > q-K/2K,q+K/a-Q> ~ ~ ~ " • (45) 

m ^n,q-K/2 — £^,q-K/2 

The two terms in the sum of Eq. (Q) are closely related to the matrix elements for radiative phonon-assisted 
exciton recombination. We show the radiative recombination processes of an exciton in Fig. 3, diagrams F and F', 
where an exciton recombines emiting a phonon of momentum Q and a photon of momentum k 7 . The similarity of 
the upper part of diagrams C, C, D, and D' with F, and F', enables us to extract information about the phonon- 
assisted Auger mechanism. Experimentally, the primary luminescence mechanism of ortho excitons in CU2O is via 
phonon-assisted recombination involving the T^ 2 optical phonon. Other phonon-assisted recombination processes of 
the ortho excitons are at least 30 times weaker. For the para excitons the only allowed optical phonon recombination 
mechanism is via a r^~ 5 phonon, and the integrated intensity is about 500 times weaker than the T^ 2 ortho-exciton 
phonon-assisted mechanism. This difference is due to the band structure of the material and the symmetry of the 
bands that are involved in the transition. Group theory allows both odd-parity bands to assist the ortho-exciton 
radiative recombination, so for the orthoexcitons both processes F and F' are allowed. In para-exciton recombination, 
in contrast, only the deep valence band can be the intermediate state (diagram F'), as illustrated in Fig. 6 p2|| . This 
analysis is valid only for non-stressed crystals; if stress is applied the mixing between the r^" and the valence bands, 
dependent on the orientation of the applied stress, is expected to affect the phonon-assisted radiative recombination 
processes. 

In the Auger problem this symmetry argument implies that for ortho-exciton recombination, the rate due to 
processes C and C is much larger than the rate due to D and D'. For para excitons, processes C and C are forbidden 
and D and D' contribute negligibly to the decay rate. To estimate the ratio of the magnitudes of the matrix elements 
for the faster ortho-exciton (Mc + Mc) and the para-exciton (Mr> + M0 1 ) Auger phonon-assisted mechanisms due 
to this effect we note that the only difference in this case as compared to the radiative recombination comes from the 
energy denominator in Eq. (E5J). An estimate of this ratio is 



M c + M C ' 



M D + M D , 



r o ,rad £c',0 — £ v,0 /t^qq £c '-° ~ £ ">° (4g) 



r p ,rad £c,0 ~ £«' ,0 £c,0 — £v' ,0 



where r 0jla d is the total decay rate of ortho excitons due to the r^~ 2 phonon and r Pjlac i is the total decay rate of para 
excitons due to the r^~ 5 phonon. As can be seen from the energy differences in the denominator of Eq. (|4^) using the 
energy levels of CU2O, \(Mc + Mc')/(Md + Md')\ ~ 3.0 x V500. The only assumption we have made in writing Eq. 
(fl6|) is the reasonable one that the oscillator strengths between the bands c', v, and c, v' are comparable. Thus, the 
relative decay rate due to processes C, C as compared with the decay rate of D, D' is 5 x 1 3 . Since, as we show 
below, all other diagrams (of the form of E in Fig. 3) for the phonon-assisted Auger mechanism can be neglected, we 
conclude that para excitons have a negligible phonon-assisted Auger decay rate (in para-para collisions), compared 
with the ortho excitons, in agreement with the observed long para-exciton lifetime when no ortho excitons are present. 
Since the T^ 2 optical phonon dominates the phonon-assisted Auger decay process of ortho excitons, we neglect the 
contribution of all other phonons. 

Expanding the overlap integral (u t , q _K/2| w c',q+K/2-Q) m powers of K — Q, using k • p perturbation theory, we see 
that the matrix element for processes C and C goes as 

I 

|K _ Q| 



M c + Mc ~ |^_ rt | 2 (0 ke -P/2-K+Q - 0k e -P/ 2 )(K - Q) ■ p v , c ,, (47) 



so the matrix element Mc + Mc can be treated as constant. We estimate it using Eqs. ( f44| ) and (45) as 
M +M W (tT.A ^ P^-(K-Q)(^,o|g r - 2 K,o) 

M c + M C ' w — (K - q) ■ — — — — > 0q -, r -, Ok e +k h +Q,K+P- (48) 

n V "P /k £oo|K-Q| 2 m (£ c / l0 - e„,o) (e c ,o - £c',o) 







The matrix element (u c i t o\H r - l M c,o) equals D r - [h/ (2pfluj r - )] 1 / 2 in magnitude, where o; r - is the zone center fre- 
quency of the specific phonon, « 13.8 meV, p is the mass density of the material, and D r - is the deformation potential 

1 12 

involving the rj~ 2 optical phonon-mediated interband transition between the T§ and the d band. Therefore, Eq. ( |48| ) 
can be written as 

1 /2 

M MM ~ -ioo- e2aB h l P ^ c 'l m€H Dr ™ I 11 I A //(q\ 

Mc + «(?« « 128?r — — — — — o ke +k h +Q,K+p- (49) 

eoofi m (£ c ',o - e«,o) [1 + (A^asr^ (e c ,o - £ c ',o) ^pfiwp- J 

The deformation potential £> r - has not been measured directly, but we can extract an estimate from Ref. p9[ . We 
use the measured para-exciton lifetime of 13 ps at a temperature of 10 K and the relative integrated intensities of 
the IY 2 ortho-exciton phonon-assisted mechanism with respect to the T^ 5 para-exciton phonon-assisted mechanism 
(approximately 500, as mentioned earlier), to fit the decay rate of the ortho excitons due to the phonon-assisted 
recombination, 26 ns at the same temperature of 10 K. This analysis implies that D T - w 2.5 eV/A. 

For the processes of the form of E shown in Fig. 3, which correspond to the diagrams B and B' (direct Auger 
process) with additional emission of a rj" 2 phonon, we have 

M E = J2 <Mk h -P/2^ff(k e q K/2) <5k e +k h +Q,K+p. (50) 

q £ C ,k h -P — £c',k e -K - IIUJQ 

Since the overlap integral is less than unity, an upper bound on Me is 



M E < 



(u c >,o\H r - \u c>0 ) 



£ c ,0 — £ c',0 q 



^2 ^qKff(k e - q) <$k e +k h +Q,K+P 



(51) 



The sum over q in Eq. ([si]), the convolution of the wave function of the relative electron- hole motion (f> Ie \(x) times 
the Coulomb interaction V e s(x), can be calculated analytically, yielding 

M < 1 A ^ ( h 4 (^) 1/2e2 x ,™ 

M ^irm^(e c , -s c ,, ) [w^r) ^ITT(fc^I^ +kh+Q - K+p - (52) 

Let us proceed to the calculation of the decay rate. If M p h is the matrix element for the phonon-assisted process, the 
decay rate of two colliding excitons with momenta K and P giving an electron of momentum k e , a hole of momentum 
kh with simultaneous emission of an optical phonon of momentum Q, is 

Tk,p = y Yl l M ph| 2 (l-™ c ,kJ(l-?VkJ<S(£K + £p -hoj Q -e cMe -e vM J. (53) 

k 6 ,k h ,Q 

We neglect the inverse process at very low lattice temperature, where no phonons are present. 

The contribution of the (dominant) processes C and C to the Auger phonon-assisted decay rate, defined in Eq. 
(|l|),is, using Eq. ©, 

r _x „, 2"xl0-' M e 4 g na s \E^l l D r;°B? 
A ' ph 7T 2h 3 el c mpa 7 B B m (e C:0 - e C ',o) 2 (e c ',o - £„,o) 2 7iw r - 

In the approximation of Eq. (47), where the matrix element is constant, T^ ph is temperature independent at low 
lattice temperature and the corresponding inverse scattering time is proportional to the density of the total number 
of excitons, T^p h ~ n. If we measure the density in units of 10 18 cm~ 3 , Eq. (|5~i| ) gives the decay rate 



T 



^wO^nns" 1 , (55) 



A,ph 

which is the order of magnitude of the observed decay rate, but not numerically accurate. 
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For process E, the approximate expression on the right side of ( p2| ) for the matrix element yields a lower bound on 
the contribution TA,ph,E to the Auger phonon-assisted scattering time, 

_ x 2 10 xl0- 4 M e 4 ^ 3 (J> r -a B ) 2 

TW ~ * 2ft 3 e*, pa% ™ B ( £c , - ec,o) a ^ r - ' ( j 

Numerically T^ ph E < 9 x 10~ 3 n ns" 1 , with n measured in units of 10 18 cm -3 , which is comparable to the decay rate 
of the para excitons, i.e., the smallest decay rate of our problem. Processes of the form of E are in fact much slower 
than this rate, since we have assumed in (|5^) that the overlap integral is unity. Furthermore, the process E', of the 
form of E with the Coulomb exchange coupling instead to the positively charged hole, enters with a change of sign, 
partially cancelling E. 



V. EFFECTS OF IDENTITY ON THE AUGER DECAY RATE 

We turn now to the problem of the identity between the electrons and holes of the two excitons. The total wave 
function of two excitons of momenta K and P with angular momentum functions \ a and \p respectively, is 

*i,tot = i[*K(l,2)xa(l,2)* P (l',2')x/3(l',2') - *K(l',2)xa(l', 2)* P (l,2')x,j(l, 2') 

-* K (1, 2')Xa(l, 2')*p(1', 2) X/3 (1', 2) + * K (1', 2')Xa(l', 2')*p(1, 2) X /j(1, 2)], (57) 

where the variables 1, V , 3 refer to electrons and 2, 2', 4 to holes. The wave fuction of the final state is 

%tot = *k e ,k h (3,4) X7 (3,4). (58) 

The calculation we have presented has been done without effects of identity taken into account. Thus, if the initial 
state in the matrix elements Ma, Ma>, Mq and Mq 1 is the first term of Eq. ([s?]), the second term of Eq. ( [57|) gives 
the initial state in the matrix elements Mb, Mb>, Md and Md>; the last two give exchange terms. Since the phonon- 
assisted processes C and C dominate the Auger collisions, we examine here the effect of identity of the electrons and 
the holes on the corresponding matrix elements Mq and Mc> ■ 

The Coulomb interaction does not produce any spin-flips, so the angular momentum functions \i factor out in all the 
matrix elements, independent of the spatial part of the calculation. For ortho-para and for ortho-ortho collisions with 
different J z , only the first or the last term of Eq. (|5^) contributes to the matrix element, because of the orthogonality 
of the angular momentum wave functions. The actual matrix element, therefore, for the case that the effect of identity 
has been taken into account is 1/2 (from the normalization constant) times the result of our calculation from the 
previous section. For ortho-ortho collisions with the same J z , both the first as well as the last term in Eq. j5^ ) 
contribute to the matrix element and therefore the numerical factor which multiplies the matrix element in this case 
is 1/2+1/2=1. 

In a gas of orthoexcitons with random J z , there are nine possible combinations, with respect to their angular 
momentum, that two ortho excitons collide. Statistically 3/9 of the collision events are between ortho excitons with 
the same J z , and 6/9 between ortho excitons with different J z . Consequently, if Ca (~ 0.5 from Eq. ([55])) is the 
proportionality factor in our calculation without taking into account the effect of identity, the effect of identity on 
the decay rate is to multiply the original rate by a factor (3/9) 1 + (6/9) (1/2) 2 = 1/2 for ortho-ortho collisions and 
(1/2) 2 for ortho-para collisions. Thus for the dominant phonon-assisted Auger process 

1 f dN °Y - 1 CA ( M 1 ^ /cq\ 

1 ^ P =^-^n 01 (60) 



N P V dt J Aph rl ph 4 
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where N a and N p are the numbers of ortho and para excitons, respectively. To determine the net loss of ortho and 
para excitons due to the Auger process, we take into account the fact that the Auger-ionized excitons reform in 
random angular momentum states, and find 

1 ' dN *V CA n Q U- n A. (62) 



N P V dt K et 16 "V n. 



VI. EXPERIMENTAL OBSERVATIONS - CONCLUSIONS 

The two basic experimental observations related to the Auger process are shown in Figs. 2 and 7. Figure 2 shows 
the total decay rate per particle of the ortho excitons as function of the ortho-exciton density. As seen from this 
graph, the total decay rate is approximately proportional (the exact power is 0.8) to the ortho-exciton density. In 
this experiment a mode-locked, cavity dumped Argon-ion laser produces nanojoule pulses with about a 100-ps length. 
For a photoluminescence time resolution of about 100 to 300 ps, this excitation pulse is effectively a delta- function 
in time, and the evolution of the system is observed without creation by the laser of further particles. The excitons 
form on a timescale shorter than the detection limits, within a few nanoseconds. The much weaker radiative efficiency 
of the para excitons makes it difficult to observe them during this short time period. As mentioned earlier, this 
experiment shows quantum saturation of the ortho excitons, i.e., they move along an adiabat (— /i /fcsT w 0.2, where 
fi a is the chemical potential of the ortho excitons) very close to the condensation line, without condensing; along 
adiabats n a ~ T 3 / 2 . Our result, Eq. (|sl|) , gives an approximate linear dependence of the Auger decay rate per ortho 
exciton as function of the ortho-exciton density. Experimentally the Auger decay rate cannot be measured directly, 
since ortho-para interconversion || contributes almost as much as the Auger to the total rate shown in Fig. 2. Since 
the decay rate due to the conversion of ortho excitons to para excitons is approximately ~ T 3 / 2 , which is ~ n Q along 
an adiabat, the contribution of the Auger decay rate to the total rate experimentally is also ~ n , in accordance 
with our theoretical calculation. In the numerical simulation in Ref. (2), the value of 0.40 ns -1 for the Auger decay 
constant Ca (with the density measured in units of 10 18 cm -3 ) reproduces the experimental results; this value should 
be compared with 0.5 ns -1 from Eq. (|55|). 

If the mode-locker is removed from the laser, the cavity-dumped mode provides 10-ns long pulses with about an 
order of magnitude more energy. Figure 7(a) shows data from lightly stressed crystals with such a long-pulse excitation 
The laser profile (triangles), the number of ortho excitons in the lowest ortho-exciton level (open circles), and 
the number of para excitons (black dots) are shown as function of time in this figure. The stress splits the triply 
degenerate ortho-exciton level into three components and only the lowest of the three is significantly populated, 
leading to a closer proximity to the condensation line than in the unstressed case. No ortho-exciton condensation is 
observed. Under these conditions the para-exciton density was determined from the relative intensities of ortho and 
para excitons, combined with the spectroscopically determined density of the ortho excitons. This analysis yields the 
striking evidence that the para excitons condense shown in Fig. 7(b), a graph of the corresponding trajectories for the 
ortho excitons (open circles) and the para excitons (black dots) in the density-temperature plane. The straight line 
here is the condensation phase boundary, which is identical for para excitons and ortho excitons in the stressed case. 
The para excitons in this case are Bose condensed for times later than 8 ns. A crucial feature of the para excitons 
in these data which is related with the Auger process is their significantly smaller decay rate than that of the ortho 
excitons. More specifically, their lifetime is on the order of 100 ns at a density 10 19 cm -3 , while the lifetime of the 
noncondensed ortho excitons (ss 0.1 nsec) is consistent with the phonon-assisted Auger and ortho-para conversion 
mechanisms discussed earlier, in the case of short-pulse excitation. The para-exciton decay time of 100 ns refers to 
late times, when the production of excitons due to the laser is negligible, the ortho-exciton density is very low, and 
essentially only para excitons exist. This observation shows that para-para Auger collisions lead to a very low decay 
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rate. As we showed in Sec. IV, the reason for this low rate is that the band with the required symmetry which assists 
the para-para phonon-assisted Auger processes is a deep valence band. See Fig. 6. In contrast, for ortho-ortho, or 
ortho-para collisions the relevant band is very close to the conduction band. In the numerical simulation in Ref. Q 
of the time dependence of excitons, we have shown that the ortho excitons move along the phase boundary without 
crossing it while the para excitons condense, as a result of the Auger heating of the othro excitons and phonon cooling 
of both the ortho and para excitons. 

The result of high uniaxial stress ~ 3.6 kbar on the Auger process has been resolved in Ref. JhJ. For comparison, 
in the data of Fig. 7 the stress is « 0.36 kbar; in the data of Fig. 2 there is no applied stress. The experiment in 
Ref. |l0| shows clearly that para-para Auger collisions are allowed in stressed crystals. In this experiment laser light 
produces the excitons away from the well, which drift a distance of « 400 /jm in the stress gradient to the shear-stress 
maximum where they are confined. The number of ortho excitons in the well should be negligible, since the relatively 
short lifetime of ortho excitons precludes travel over 400 /im distances; also the expected number of ortho excitons 
due to thermal excitations is negligibly small under the specific experimental conditions (T w 2 K). Experimentally, 
however, the measured number of ortho excitons is approximately proportional to the square of the para-exciton 
number, and the maximum para-exciton density goes as the square root of the laser power. These observations have 
been attributed to the Auger decay of para excitons . In this case that the crystal is highly-stressed, the symmetry 
selection rules we have used in both the phonon-assisted, as well as the direct processes are expected to break down. 
An example of this symmetry breaking is that the (zero stress-forbidden) direct recombination line of para excitons 
has been observed clearly in stressed crystals; as far as we know, however, the result of stress on phonon-assisted 
recombination has not been investigated yet. 

Stress strongly influences the Auger process, starting with the phonon-assisted mechanism; both the dominant 
phonon-assisted Auger process involving the r^ 2 phonon as well as the direct Auger become allowed for para-para 
collisions in stressed crystals, for specific directions of the stress. More specifically, group theory predicts [^0) that if 
the uniaxial stress is along the C2 or the C3 crystallographic axes, the mixing of the r^" with the band allows the 
IY 2 LO-phonon to assist the process. In the specific experiment of Ref. the applied stress is primarily along the 
C4 axis, but it also has a component in the C2 direction. Moreover, uniaxial stress makes the direct Auger process 
(dominated by the processes shown in diagrams A and A') weakly allowed. Both processes are responsible for the 
presence of ortho excitons in the quantum well, which are generated by the reformation of Auger-ionized para excitons 
in random angular momentum states. In both cases the orientation of the uniaxial stress is very important, since it 
determines the selection rules in the deformed crystal Q • 

In summary, we have calculated the direct and phonon-assisted Auger decay rate, incorporating the band structure 
of CU2O, which plays an essential role because of the following effects: 1) the same (positive) parity of the two bands, 
2) the fact that the valence band is not a pure spin state and 3) the symmetry and the location of the negative-parity 
bands with respect to the conduction and valence bands which form the yellow excitons. The inverse Auger decay 
times for the direct as well as the much faster phonon-assisted Auger mechanism are proportional to the exciton 
density. Finally, we have shown that the para-exciton recombination vertex is either forbidden (direct processes), or 
negligible (phonon-assisted processes) , and consequently para-para Auger collisions have a negligible Auger decay rate 
in non-stressed crystals. 
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FIG. 1. Two classes of diagrams for the direct Auger non-radiative recombination. Time progresses from left to 
right. The initial state contains two excitons of momenta K and P, and the final state contains an ionized electron 
and hole with momenta k e and kh, respectively. The dashed line denotes the Coulomb interaction. Processes A and 
A' are the dominant ones and they are nonzero only if the recombined particle is an ortho exciton. 

FIG. 2. The total decay rate @] of ortho excitons versus the ortho-exciton density in unstressed crystals. Figure 
reproduced with the kind permission of J.P. Wolfe. 

FIG. 3. The phonon-assisted Auger process (diagrams C, C', D, D' and E), and the phonon-assisted radiative 
recombination process (F and F'). The phonon is denoted by a dashed line. Processes D and D' are the dominant 
ones because of the band structure of CU2O. In these, the recombined exciton can only be an ortho exciton. Similarly, 
diagram F is dominant compared to F'; the recombined exciton can only be an ortho exciton. 

FIG. 4. Virtual annihilation diagram of an exciton, possible only for pure spin-singlet excitons. 

FIG. 5. Schematic band structure of CU2O showing the conduction Tq band and the r^, valence bands (split 
by the spin-orbit splitting) which form the yellow and green exciton series, respectively. 
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FIG. 6. Phonon-assisted radiative recombination processes. The ortho excitons are allowed to participate in both, 
but for the para excitons only the one on the right side is allowed. The bands c, v are the rjj" and the bands which 
form the yellow excitons and are separated by 2.17 eV. The bands d and v' are the only odd-parity bands in CU2O. 
The energy difference between d and c is « 449 meV and between v and v' is rs 5.6 eV. The proximity of band d to 
band c makes both the radiative recombination p2| and the phonon-assisted Auger decay processes much faster for 
ortho excitons than para excitons. 

FIG. 7. (a) Data from lightly stressed crystals with long-pulse (10-ns) excitation ||. The laser profile (triangles), 
the number of ortho excitons in the lowest ortho-exciton level (open circles), and the number of para excitons (black 
dots) as function of time. The para excitons show a significantly smaller decay rate, (b) Corresponding trajectories for 
ortho excitons (open circles) and para excitons (black dots) in the density-temperature plane. The straight line is the 
condensation phase boundary, which is identical for para excitons and ortho excitons in the stressed case. Note that 
the para excitons are in the condensed region at times later than 8 ns. Figure reproduced with the kind permission 
of J.P. Wolfe. 
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